AD-A212  290 


DTff!  Fi i  p  cnp* 


a  j?u.  3'^-/-" 

& 


Dynamics  of  Separation  and  Reattachment  in  a  Mach  5 
Compression  Ramp-Induced  Shock  Wave  Turbulent 
Boundary  Layer  Interaction 


Final  Report 


R.  A.  Gramann  and  D.  S.  Dolling 


July  31,  1989 

U.  S.  Army  Research  Office 

Grant  DAAL03-86-G-0045 

The  University  of  Texas  at  Austin 
Aerospace  Engineering  and  Engineering  Mechanics  Dept. 


Approved  for  Public  Release; 
Distribution  Unlimited. 


'  jL 


1».  REPORT  SECURITY  CLASSIFICATION 


iinrafT-T-iurr.i 


2a.  SECURITY  CLASSIFICATION  AUTHORITY 


2b.  DECLASSIFICATION/ DOWNGRADING  SCHEDUL 


4.  PERFORMING  ORGANIZATION  REPORT  NUMBER(S) 


REPORT  DOCUMENTATION  PAGE 


6a.  NAME  OF  PERFORMING  ORGANIZATION 
University  of  Texas  at  Austin 


6c  ADDRESS  (Gty,  Start,  and  ZIP  Cod*) 

Aerospace  Engineering  Department 
WFW  217,  Austin,  TX  78712 


6b.  OFFICE  SYMBOL 
(If  applicable) 


8a.  NAME  OF  FUNDING /SPONSORING 
ORGANIZATION 

U.  S.  Army  Research  Office 


8c  ADDRESS  (City.  Start,  and  ZIP  Cod*) 

P.  0.  Box  12211 
Research  Triangle  Park,  NC 


8b.  OFFICE  SYMBOL 
(If  applicable) 


27709-2211 


lb.  RESTRICTIVE  MARKING 


3  .  DISTRIBUTION /AVAILABILITY  OF  REPORT  s 

Approved  for  public  release; 
distribution  unlimited. 


S.  MONITORING  ORGANIZATION  REPORT  NUMBER(S) 

0  J37L  Z.l-ei-P 


7a.  NAME  OF  MONITORING  ORGANIZATION 
U.  S.  Army  Research  Office 


7b.  ADDRESS  (Gty,  State,  and  ZIP  Cod*) 

P.  0.  Box  12211 

Research  Triangle  Park,  NC  27709-2211 


9.  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 


10.  SOURCE  OF  FUNDING  NUMBERS 


PROGRAM 
ELEMENT  NO. 


PROJECT 

NO. 


WORK  UNIT 
ACCESSION  NO. 


11.  TITLE  (Include  Security  Classification) 

DYNAMICS  OF  SEPARATION  AND  KEATTAChMEOT  IN  A  MACH  5  COMPRESSION  RAMP  ET£W. 


12.  PERSONAL  AUTHOR(S) 

R.  Granann,  D.  S.  Dolling 


13a.  TYPE  OF  REPORT  Il3b.  TIME  COVERED 

Final  I  from  6/1/86  to 


14.  DATE  OF  REPORT  (Year,  Month,  Day) 

August  7,  1989 


.  PAGE  COUNT 

5 


16.  SUPPLEMENTARY  NOTATION  _ 

The  view,  opinions  and/or  findings  contained  in  this  report  are  those 
of^Jhe  authgr(|) .and  should  not jbe4 const jugd  as  an  official  Department  of  the  Army  position. 


COSATI  CODES 


GROUP  SUB-GROUP 


18.  SUBJECT  TERMS  ( Continue  on  reverie  if  necessary  and  identify  by  block  number) 

Shock  Wave  Turbulent  Boundary  Layer  Interaction 
Unsteady  Flow 
Turbulent  Separation 


19.  ABSTRACT  ( Continue  on  reverse  If  necessary  and  identify  by  block  number) 

The  separation  process  in  cylinder vand  unswept  compression  ramp  induced  shock  wave  turbulent 
boundary  layer  interactions  -has  been^examined  using  a  conditional-analysis  of  wall  pressure 
fluctuation  measurements.  The  tests  were  performed  in  a  Mach  5  blowdown  tunnel  under 
adiabatic  wall  terrperature  conditions .  The  conditional  analysis^lrasishown  that  the 
instantaneous  separation  position  is  at  or  close  to  the  instantaneous  separation  shock  foot 
in  these  interactions.  The  separation  line  as  indicated  by  traditional  surface  tracer 
methods  is  at,  or  close  to,  the  downstream  boundary  of  a  region  of  intermittent  separation. 
The  dynamics  of  the  separation  bubble  in  the  unswept  oonpression  ramp  flcwfieldfihave'also 
-beert>examined  and  initial  results  suggest  that  the  motion  of  the  separation  point  and 
reattachment  point  locations  are  correlated.  The  preliminary  results  indicate  the 
reattachment  point  is  at  its  downstream  locations  when,  the  separation  shock  is  at  its 
upstream  locations.  I]  r  ,  .  r  -  L,  -4  S'sn  •  Sir 


20.  DISTRIBUTION /AVAILABIUTY  OF  ABSTRACT 
□  UNCLASSIFIED/UNUMITED  □  SAME  AS  RPT. 


22a.  NAME  OF  RESPONSIBLE  INDIVIDUAL 


□  DTIC  USERS 


*  ’  «*—  *  t  '  i  r  /  r  r 

:  /  / ;/  r  bit -l'  *  •  f .  /*  .*  *  -  * 

/">  /  J  *  i  Jo*'  .  *  . 

:  .nr*  7  ?  rscl  -  .■■■* v.':  ■.  : 


21.  ABSTRACT  SECURITY  CLASSIFICATION 
Unclassified 


22b.  TELEPHONE  (Include  Area  Code)  22c.  OFFICE  SYMBOL 


DO  FORM  1473, 84  mar 


83  APR  edition  may  be  used  until  exhausted. 
All  other  editions  are  obsolete. 


SECURITY  CLASSIFICATION  OF  THIS  PAGE 
UNCLASSIFIED 


TABLE  OF  CONTENTS 


Section 

•4*. 

Page 

List  of  Figures 

ii 

List  of  Tables 

iii 

List  of  Appendices 

iii 

Statement  of  Problem 

1 

Summary  of  Results 

2 

List  of  all  Publications 

24 

List  of  Personnel 

24 

References 

25 

Appendices 

26 

LIST  OF  FIGURES 


Figure  Title  Page 

1  Sideview  of  Flowfield  with  Instrumentation  12 

2  Separated  Cross  Correlations, Cylinder  and  Ramp  13 

3  Conditional  Cross  Correlations,  Circular  Cylinder  14 

4  Conditional  Cross  Correlations,  Compression  Ramp  15 

5  Separated  Cross  Correlations,  Compression  Ramp  16 

6  Time- Averaged  Pressure  Distribution  17 

7  Wall  Pressure  RMS  Distribution  18 

8  Cross  Correlations  Upstream  of  ’R’  19 

9  Cross  Correlations  Near  ’R’  20 

10  Cross  Correlations  Downstream  of  ’R’  21 

1 1  Cross  Correlations  at  Mid-Ramp  Location  22 

12  Simplified  Flowfield  Model  23 

13  Shear  Stress  Model  35 


n 


LIST  OF  TABLES 


Table  Title  Page 

•u. 

1  Freestream  and  Boundary  Layer  Conditions  1 1 


LIST  OF  APPENDICES 


Apppendix  Title 

A  AIAA  Paper  88-4676,  AIAA  Measurements  and  Instrumentation 

Conference,  Atlanta,  GA,  September  1988,  "Detection  of  Shock 
Induced  Separation  using  Wall  Pressure  Fluctuations" 

B  Surface  Tracer  Response  Explanation-from  AIAA  Paper  86-1033, 

AIAA  Fluids,  Plasma  Dynamics,  and  Lasers  Conference, 

Atlanta  GA,  May  1986. 


m 


1 


Statement  of  Problem 

The  purpose  of  this  study  was  to  investigate  the  dynamics  of  the  separation 
bubble  in  an  unswept  compression  ramp-induced  shock  wave  turbulent 
boundary  layer  interaction.  This  flowfield  has  been  examined  numerous 
times  in  earlier  studies,  with  a  variety  of  techniques[l-7].  These  studies  have 
addressed  several  features  of  this  flowfield,  such  as  spanwise  variations  (three 
dimensional  structures/effects),  the  separation  shock  motion,  and  mean 
features  of  the  flowfield,  including  the  overall  structure.  However,  little  time 
dependent  information  on  the  separation  process,  the  separation  bubble,  or 
the  region  of  reattachment  has  been  obtained.  To  develop  a  complete 
understanding  of  the  interaction,  information  on  these  latter  aspects  of  the 
flow  must  be  known.  In  particular,  information  on  the  separation  bubble  is 
needed  to  complete  the  overall  flowfield  picture.  With  a  complete  time- 
dependent  description  of  the  interaction,  insight  into  the  physical 
mechanisms  controlling  this  interaction  will  be  available.  This  insight  will 
also  provide  guidance  for  the  required  improvement  in  the  numerical  models 
being  used  to  predict  these  interactions. 


2 


Summary  of  Results 

To  date,  several  interesting  results  have  been  obtained  through  the  test  and 
analysis  program  performed  under  this  grant.  Some  of  these  have  already 
been  published  and  consequently  will  only  be  summarized  here.  Continuing 
work,  needed  for  answering  one  or  two  remaining  questions  will  also  be 
described. 

First,  a  method  of  determining  flow  direction  using  wall  pressure  fluctuations 
has  been  developed  and  has  been  used  to  detect  separated  flow  in  two 
interactive  flowfields.  Two  channels  of  wall  pressure  fluctuations  from  two 
closely  spaced  high  frequency  pressure  transducers  are  used  in  the  method. 
Thus  the  method  is  non-intrusive  and  relatively  easy  to  use.  Kulite  miniature 
pressure  transducers  are  mounted  flush  with  the  test  surface,  in  the  region  of 
the  interaction  of  interest.  Figure  1  shows  a  centerline  view  of  one 
interaction  and  one  typical  instrumentation  arrangement.  In  this  figure,  the 
transducers  are  located  in  the  intermittent  region,  which  is  characterized  by 
large  pressure  fluctuations  due  to  the  separation  shock  oscillations  within  this 
region.  The  percentage  time  that  a  position  in  the  intermittent  region  is 
downstream  of  the  separation  shock  is  known  as  that  location’s 
"intermittency"  (y)  value.  The  upstream  boundary  of  the  intermittent  region  is 
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the  upstream  boundary  of  the  interaction  (y=0%).  The  downstream  boundary 
of  the  intermittent  region  (Y=100%)  is  typically  delineated  by  the  separation 
line  determined  using  surfacejracers.  Downstream  of  this  point,  the  flow  is 
fully  separated  at  all  times. 

Continuous  data,  sampled  simultaneously  on  the  two  channels  at  rates 
ranging  from  200  to  500  kHz  are  stored  in  CPU  memory  and  then  transferred 
to  disk.  Analysis  is  then  performed  on  the  data.  By  cross  correlating  selected 
data  corresponding  to  a  certain  flow  conditions,  features  of  that  flow,  such  as 
flow  direction  and  time-averaged  large  scale  structure  velocity,  can  be 
determined.  This  analysis  method,  described  in  detail  in  Publication  2(see  p. 
2 U),  was  used  to  determine  the  features  of  the  separated  region  in  the  28° 
compression  ramp  induced  interaction.  The  freestream  Mach  number  was 
4.95.  Wall  conditions  were  nearly  adiabatic;  the  boundary  layer  on  the  tunnel 
floor  developed  naturally  and  was  fully  turbulent  at  the  test  locations. 
Additional  flow  conditions  are  given  in  Table  1. 

Separation  process: 


The  analysis  method  described  briefly  above  has  been  used  to  examine  the 
separation  process  in  a  shock  wave  turbulent  boundary  layer  interaction. 
Data  from  a  cylinder  induced  interaction  and  a  compression  ramp  interaction 
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have  been  analyzed  and  several  interesting  points  concerning  separation  have 
been  clarified.  Figure  2  shows  cross  correlation  results  obtained  using  the 
signals  from  two  pressure  transducers  placed  just  downstream  of  ’S,’  the 
separation  line  determined  using  surface  tracers.  The  transducers  were 
mounted  in  the  streamwise  direction,  spaced  0.292  cm  apart.  Both  the 
cylinder  and  ramp  cross  correlations  show  a  distinctive  double  peak, 
indicative  of  separated  flow.  A  more  detailed  discussion  of  this  result  is 
presented  in  Appendix  A. 

To  determine  flow  direction  in  the  intermittent  region,  data  corresponding  to 
flow  downstream  of  the  separation  shock  were  extracted  from  continuous 
time  data  obtained  with  both  transducers  mounted  in  the  intermittent  region. 
Upstream  of  the  separation  shock,  undisturbed  boundary  layer  flow  exists. 
Hence  die  direction  anu  broadband  velocities  for  this  portion  of  the  flow  are 
known.  The  "Conditionally  Extracted  Analysis'  Data  Sets"  (CEADS) 
corresponding  to  flow  downstream  of  the  separation  shock  were  analyzed 
using  the  cross  correlation  algorithm.  Using  the  cross  correlations,  the  flow 
direction  and  velocities  downstream  of  the  shock  can  be  calculated.  Figure  3 
shows  several  results  from  locations  in  the  intermittent  region  of  a  cylinder 
induced  interaction.  The  "double  peak"  pattern,  indicative  of  separated  flow, 
is  clear  and  has  the  same  timing  as  those  in  Figure  2.  Thus  separation  occurs 
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immediately  downstream  of  the  separation  shock  at  all  locations  within  the 
intermittent  region  of  the  cylinder  induced  interaction. 

Similar  analyses  of  ramp  induced  interaction  data  were  performed  as  well. 
Typical  results  are  shown  in  Fig  4.  Again,  separated  flow  characteristics  are 
seen  in  the  cross  correlation  at  several  locations  within  the  intermittent 
region.  Thus  separation  occurs  within  the  intermittent  region,  immediately 
downstream  of  the  oscillating  separation  shock.  Both  cylinder  and  unswept 
compression  ramp  induced  interactions  exhibit  this  feature.  Therefore,  the 
intermittent  region  is  a  region  of  intermittent  separation.  Further,  the 
separation  line.’S,’  is  close  to  the  downstream  boundary  of  this  region  of 
intermittent  separation.  A  physical  explanation  of  this  phenomenon  is 
provided  in  Appendix  B. 

Separation  Bubble  Dynamics: 

Data  have  been  obtained  in  the  separated  flow  (i.e.  downstream  of  ’S’)  of  the 
unswept  compression  ramp  induced  interaction,  both  upstream  of  the  ramp 
comer  and  on  the  ramp  face.  The  data  just  upstream  of  the  ramp  comer  have 
been  used  to  characterize  the  separated  flow  on  the  ramp  face.  Figure  5 
shows  a  cross  correlation  result  obtained  with  both  transducers  closer  to  the 
ramp,  along  with  the  previous  result  shown  in  Fig.  2.  This  result  is 
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characteristic  of  separated  flow  near  the  ramp  comer  and  will  be  used  for 
comparison  with  cross  correlations  on  the  ramp  face. 

Data  near  reattachment  have  been  obtained  for  two  separation  shock 
positions,  as  well  as  continuous  time.  800  records  total  are  obtained  during 
the  continuous  time  data  acquisition,  providing  approximately  one  second  of 
continuous  time  data.  These  data  have  been  analyzed  to  determine  if 
information  regarding  the  the  separated  flow  can  be  obtained  near 
reattachment.  Based  on  data  obtained  for  fixed  separation  shock  position, 
the  preliminary  results  indicate  that  reattachment  (or  the  flow  character  near 
reattachment)  is  dependent  on  the  shock  position.  Figures  6  and  7  show  the 
mean  pressure  and  RMS  of  the  wall  pressure  fluctuations  distributions  within 
the  interaction.  The  time-averaged  results  are  similar  to  results  obtained  by 
other  investigators  in  other  facilities[l-5J.  It  is  interesting  to  note  that  the 
mean  pressure  and  RMS  for  the  furthest  upstream  separation  shock  location 
(y  =  0-8%)  are  significantly  less  than  the  continuous  time  values.  When  the 
separation  shock  is  between  y  =  30%  and  y  =  40%,  the  mean  pressure  and 
RMS  are  approximately  equal  to  the  continuous  time  pressures.  This  suggests 
that  the  wall  pressure  and  its  fluctuations  on  the  ramp  face  are  a  function  of 
separation  shock  position.  Further  data  are  currently  being  obtained  to 
expand  the  regions  measured  on  the  ramp  face  and  increase  the  number  of 


7 


shock  positions  examined. 

Cross  correlations  were  calculated  from  data  from  three  pairs  of  transducers 
located  near  ’R,’  the  reattachment  line  determined  using  surface  tracers. 
Results  at  one  mid- ramp  location,  appoximately  1.6  8q  downstream  of  ’R,’ 
were  also  obtained.  Figures  8-11  show  the  resulting  cross  correlations  of  the 
wall  pressure  fluctuations  for  each  position.  In  each  figure,  both  continuous 
time  data  and  specific  shock  position  data  (or  "shock-fixed")  data  are 
presented. 

Figure  8  shows  the  results  with  both  transducers  located  upstream  of  the 
reattachment  line,  ’R.’  Although  a  change  in  magnitude  in  the  correlation 
coefficients  is  seen,  the  character  of  the  cross  correlations  from  the  three  data 
sets  are  the  same.  Figure  9  shows  cross  correlations  of  data  from  two 
transducers  straddling  ’R.’  A  change  in  character  is  seen  between  the 
continuous  data  result  and  the  "shock-fixed"  results.  Specifically,  the  distinct 
dips  on  either  side  of  the  dominant  peak  are  present  in  the  "shock-fixed"  data, 
but  are  barely  perceptable  at  -x  for  continuous  time  data.  Figure  10  shows  the 
same  trend.  The  dips  are  less  pronounced  for  the  "shock-fixed"  data,  and  no 
dip  is  seen  in  the  continuous  time  result.  Both  transducers  in  this  case  are 
mounted  close  to,  but  downstream  of  ’R.’  The  "dips"  are  indicative  of 
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separated  flow.  Further  downstream  on  the  ramp  face,  separated  flow  exists 
for  a  smaller  fraction  of  the  time.  Thus  the  continuous  time  signals  show  less 
and  less  of  the  separated  characteristics  (i.e.  cross  correlation  "dips")  as 
distance  downstream  of  the  ramp  comer  is  increased.  However,  the  shock- 
fixed  data  clearly  show  the  separated  character  exists  at  positions 
downstream  of  ’R.’  At  these  locations,  little  indication  of  separated  flow  in 
the  continuous  data  cross  correlation  exists. 

Figure  11  shows  the  correlations  of  signals  from  transducers  mounted  at 
mid-ramp  locations.  No  change  in  the  cross  correlations  for  the  3  cases  (i.e. 
continuous  time,  y  =  0-8%,  and  y  =  30-40%)  is  seen.  This  confirms  the 
sensitivity  of  the  method  for  the  locations  closer  to  the  ramp  comer.  These 
data  suggest  the  reattachment  point  motion  is  correlated  with  the  separation 
shock  motion.  If  no  correlation  with  shock  position  existed,  the  shock-fixed 
data  would  repeat  the  continuous  time  data,  since  random  data  sets  (i.e.  data 
sets  corresponding  to  different  flow  conditions)  would  be  used  to  used  to 
calculate  the  shock-fixed  cross  correlations.  The  continuous  time  data  has  all 
flow  conditions  for  a  given  position  represented  in  its  results.  Since  a  change 
occurs  when  shock-fixed  data  is  used  only  certain  flow  conditions  are 
occuring  for  a  given  separation  shock  location.  Additional  data  are  needed  to 
expand  the  information  available,  both  for  more  shock  positions  and  ramp 
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locations.  These  data  necessary  to  verify  these  conclusions  are  currently 
being  acquired. 

Thus,  preliminary  results  suggest  the  flow  character  at  locations  near 
reattachment  are  dependent  on  separation  shock  position  (or  separation  point 
position).  The  results  suggest  that  when  the  separation  shock  is  upstream,  the 
reattachment  point  is  downstream. 

The  shock  motion/reattachment  point  motion  result  is  further  substantiated 
by  turning  angle  calculations  in  the  intermittent  region  based  on  the  pressure 
rise  across  the  separation  shock.  The  pressure  rise  turning  angles  indicate, 
using  a  straight  line  approximation  for  the  separation  bubble  boundary,  that 
the  reattachment  point  moves  approximately  0.1  inch  with  separation  shock 
motion  ranging  from  y  =  0 %  to  y  =  100%.  The  reattachment  point  is  at  its 
upstream  location  when  the  shock  is  downstream  (high  y),  and  vice  versa. 
Figure  12  shows  a  sketch  of  the  approximated  flowfield. 

Summary: 


The  separation  process  in  cylinder  and  unswept  compression  ramp  induced 
shock  wave  turbulent  boundary  layer  interactions  has  been  examined. 
Separation  occurs  across  the  separation  shock  in  these  interactions.  Thus  the 
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intermittent  region  is  also  a  region  of  intermittent  separation.  The  separation 
line  determined  using  surface  tracer  methods  is  at,  or  very  close  to,  the 
downstream  boundary  of  this  intermittent  separation  region. 

The  separation  bubble  in  an  unswept  compression  ramp  flowfield  has  been 
examined  and  initial  results  suggest  the  separation  point  and  reattachment 
point  locations  are  correlated.  Cross  correlations  near  reattachment  are 
dependent  on  separation  shock  position.  The  preliminary  results  indicate  the 
reattachment  point  is  at  its  downstream  locations  when  the  separation  shock 
is  at  its  upstream  locations.  Further  information  on  the  downstream  shock 
locations  will  be  obtained  shortly,  completing  the  data  base. 

The  time  dependent  pressure  values  in  the  intermittent  region  and  on  the 
ramp  face  are  substantially  different  from  the  time-averaged  values.  This 
suggests  that  the  computational  models  being  used  to  calculate  these 
flowfields  must  incorporate  these  time  dependent  phenomena  before  adequate 
results  will  be  obtained.  The  mean  results  are  not  an  sufficiently  accurate 


representation  of  the  flowfield. 


Wind  Tunnel  Flow  Conditions 


Parameter 

Tunnel  Floor 

M 

oo 

4.90  ±.02 

Uoo 

741  m/s  (2432  ft/s) 

Re 

OO 

53.3  x  106  m'1  (16.2  x  10^  ft*1) 

To 

330K  (595‘R) 

Po 

2.09  x  106  N/m2  (304  psi) 

X 

0.74  m  (29  in)  from  throat 

5o 

1.62  xl0'2m  (0.63  in) 

8* 

5.23  x  10’3  m  (0.206  in) 

e 

4.54  xl04m  (1.83  xlO"2  in) 

n 

0.115 

Ree 

23.4  x  103 

mm 

9.9  x  104 

Table  1  Freestream  and  Boundary  Layer  Conditions 
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Figure  1  -  Sideview  of  Flowfield  with  Instrumentation  Installed 
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Cross  Correlation  of  Channels  1  &  2 
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Figure  2  -  Separated  Cross  Correlations 
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Figure  3  -  Conditional  Cross  Correlations,  Circular  Cylinder 
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Figure  4  -  Conditional  Cross  Correlation,  Compression  Ramp 
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Figure  5  -  Separated  Cross  Correlations-Compression  Ramp 
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Figure  6  -  Time-Averaged  Pressure  Distribution 
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Figure  7  -  Wall  Pressure  RMS  Distribution 
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Cross  Correlation  of  Channels  3  &  4 
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Analysis  of  Reattachment  Locations 
Upstream  Ramp  Transduce*  at  X  =  0.335  inches  (#3) 
Downstream  Ramp  Transducer  at  X  =  0.45  inches  (#4) 
Sampling  frequency  -  100000.0 


Figure  9  -  Cross  Correlations  Near  ’R’ 
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Analysis  of  Reattachment  Locations 


Upstream  Ramp  Transducer  at  X  =  0.425  inches  (#3) 
Downstream  Ramp  Transducer  at  X  =  0.54  inches  (#4) 
Sampling  frequency  -  100000.0 


Figure  10  -  Cross  Correlations  Downstream  of  ’R’ 
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Upstream  Ramp  Transducer  at  X  =  1.57  inches  (#1 1) 
Downstream  Ramp  Transducer  at  X  =  1.82  inches  (#12) 
Sampling  frequency  -  100000.0 


Figure  1 1  -  Cross  Correlations  at  Mid-Ramp  Location 


Figure  12  -  Simplified  Flowfield 
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Abstract 

A  technique  for  detecting  intermittent  shock-inducecfturbulent 
boundary  layer  separation  has  been  developed  and  tested  in  a 
Mach  5  blowdown  tunnel.  The  interaction  was  generated  by 
"semi-infinite"  circular  cylinders.  The  method  employs  two 
miniature  pressure  transducers  oriented  streamwise  and 
installed  flush  with  the  test  surface.  Through  cross  correlations 
of  the  conditionally  sampled  signals  of  the  two  transducers 
under  the  moving  shock  it  has  been  shown  that  the  flow 
downstream  of  the  instantaneous  shock  position  is  separated. 
The  results  indicate  that  in  these  flows,  the  separation  location 
indicated  by  surface  tracers,  such  as  the  kerosene  lampblack 
method,  is  actually  the  downstream  boundary  of  a  region  of 
intermittent  separation. 

Introduction 

Surface  tracer  techniques,  such  as  the  kerosene 
lampblack  method,  are  widely  used  in  high  speed  flows  to 
find  "separation  lines"  or  "lines  of  coalescence," 
particularly  in  shock  wave  turbulent  boundary  layer 
interactions^].  These  methods  are  relatively  easy  to  use 
and  produce  highly  defined,  repeatable  "separation  lines".  In 
the  case  of  the  kerosene-lampblack  method,  in  which  the 
pattern  is  lifted  off  the  surface  on  large  sheets  of  transparant 
tape,  full  scale  undistorted  records  are  obtained. 
Measurements  of  angles  and  length  scales  are  easily  made 
from  these  patterns  and  are  widely  used  for  comparison  with 
numerical  simulation  results. 

In  many  shock  wave  boundary  layer  interactions,  wall 
pressure  fluctuation  measurements  have  shown  that  the 
separation  shock  is  unsteady,  generating  an  intermittent  wall 
pressure  signal[2-7J.  A  typical  example,  in  a  Mach  3  blunt  fin 
interaction,  is  shown  in  Figure  1.  This  region  of  shock  motion 
is  known  as  the  intermittent  region.  Intermittency,  y,  is 
defined  as  the  fraction  of  time  a  pressure  transducer  is 
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downstream  of  the  shock,  and  is  an  indication  of  location 
within  the  intermittent  region.  The  intermittent  region  extends 
from  where  the  incoming  flow  is  first  disturbed  by  the  shock, 
to  close  to  the  separation  line,  ’S’,  indicated  by  surface 
tracers. 


M  =  3 


alupstream  boundary  layer 


i  f** .  *'  *  ‘  i*  *  . . .  r«  >  »  j  i  wT~  w~ 

oVoO  o'.  04  o'.  OB  o'.  1  ?  0. 16 

t  (SECONDS)  xicr1 

Figure  1  Kerosene  Lampblack  Flow  Visualization 
and  Wall  Pressure  Fluctuations 
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'In  such  flows  where  the  separation  shock  is  highly 
unsteady,  the  physical  meaning  of  these  surface  tracer  lines 
has  recently  come  into  question.  In  particular,  does  backflow 
actually  occur  upstream  of  the  surface  tracer  line?  This  might 
occur  since  the  surface  tracer  material  responds  to  the  mean 
wall  shear  stress,  and  the  technique  has  essentially  zero 
frequency  response.  Further,  the  mean  wall  shear  stress  at  a 
point  is  the  result  of  two  flow  fields  which  are  present  in  the 
intermittent  region  (i.e.  the  undisturbed  flow  upstream  of  the 
shock,  and  the  "disturbed"  flow  downstream  of  the  shock). 
What  is  needed  to  understand  what  the  surface  tracer  lines  in 
this  region  actually  represent  are  instantaneous  flow  direction 
measurements  close  to  the  surface.  Unlike  incompressible 
flow,  where  instantaneous  flow  direction  measurement 
techniques,  such  as  thermal  tufts,  are  reasonably  well 
developed[8,9],  no  relatively  straightforward,  measurement 
techniques  have  been  developed  for  high  speed  flows. 
Therefore,  the  need  for  a  relatively  simple  method  of 
detecting  "instantaneous"  flow  direction  is  clearly  evident. 

The  objective  of  the  work  reported  in  this  paper  was  to 
determine  if  flow  direction  could  be  deduced  from  wall 
pressure  fluctuations.  A  method  of  doing  this,  using  high 
frequency  response  pressure  transducers,  and  standard 
signal  conditioning  instrumentation  has  been  developed  and 
tested  in  a  Mach  5  shock  wave  turbulent  boundary  layer 
interaction.  Although  considerable  care  is  needed  in 
transducer  installation,  calibration  and  use,  such 
measurements  are  non-intrusive,  and  can  be  made  relatively 
easily  and  routinely  in  high  speed  flows.  The  equipi..wiii, 
technique,  analysis  involved,  and  some  results  are  presented  in 
this  paper. 

Experimental  Program 
Wind  Tunnel  and  Test  Conditions 

All  data  were  obtained  on  the  tunnel  floor  of  the 
University  of  Texas  Blowdown  Wind  Tunnel  under  essentially 
adiabatic  wall  temperature  conditions.  The  facility  has  a  17  x 
15  cm  test  section  and  operates  at  a  nominal  freestream  Mach 
number  of  4.9.  The  boundary  layer  developed  naturally  and 
was  fully  turbulent  at  the  test  location.  Table  1  gives  the 
incoming  boundary  layer  and  freestream  properties  as 
deduced  from  pitot  surveys.  The  models  used  for  the  study 
were  circular  cylinders,  1.27  and  1.9  cm  in  diameter,  8.9  cm 
and  7.6  cm  high  respectively.  Based  on  the  criterion  of  Ref. 
10,  both  cylinders  were  effectively  semi-infinite.  The  position 
of  the  cyl:  iders  could  be  varied  relative  to  the  fixed  location 
of  the  instrumentation  plug  described  below  so  that  different 
regions  of  the  flow  field  could  be  examined  (Fig.  2). 


Wind  Tunnel  Flow  Conditions 

Parameter 

Tunnel  Floor 

M« 

4.90  ±.02 

741  m/s  (2432  ft/s) 

53.3  xl06m'’ (16.2  xIO6  ft'1) 

To 

330' K  (595'R) 

Po 

2.09  xl06N/m2  (304  psi) 

X 

0.74  m  (29  in)  from  throat 

5o 

1.62  x  10‘2m  (0.63  in) 

8* 

5.23  xl0'3m  (0.206  in) 

8 

4.54  x  10^  m  (1.83  x  10'2  in) 

n 

0.115 

Re0 

23.4  x  103 

Cf 

9.9  x  10"4 

Table  1  Freestream  and  Boundary  Layer  Conditions 

KULTTE  TRANSDUCERS 


Figure  2  Model  and  Coordinate  System 
Instrumentation 

A  circular  instrumentation  plug  was  installed  flush 
with  the  tunnel  floor  approximately  0.74  m  from  the  nozzle 
throat.  The  cylinders  were  mounted  a  short  distance 
downstream  of  the  instrumentation  plug  location.  Wall 
pressure  measurements  were  made  using  miniature  high 
frequency  pressure  transducers.  Kulite  models  XCQ-062-15A 
or  XCQ-062-50A  transducers  were  used  for  all  tests.  These 
models  have  a  full  scale  range  of  15  psia  and  50  psia  with 
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nominal  sensitivities  of  13  mV/psi  and  2  mV/psi  respectively. 
Both  models  have  a  pressure  sensitive  diaphragm  0.071  cm 
(0.062  in)  in  diameter,  with  a  fully  active  Wheatstone  bridge 
bonded  to  it  The  estimated  frequency  response  of  the 
transducers,  with  covers  protecting  th6  diaphragm  installed, 
is  approximately  50  kHz.  Calibration  of  the  transducers  was 
performed  statically.  Earlier  work  [11]  has  shown  that  static 
calibrations  are  within  a  few  percent  of  dynamic  calibrations. 
In  all  cases,  the  transducers  were  mounted  flush  with  the  plug 
surface  with  a  streamwise  spacing  of  0.292  cm  center-to- 
center. 

Signals  from  the  pressure  transducers  were 
amplified  with  a  gain  of  200-500  and  analog  filtered  at  50 
kHz.  70-400  records  per  channel  of  data  (1  record  =  1024  data 
points)  were  then  digitized  by  a  12  bit  A/D  converter  (0-10 
volts  input)  at  sampling  rates  of  200,  250,  333,  and  500  kHz 
per  channel  and  stored  on  magnetic  tape.  All  data  acquisition 
and  subsequent  analysis  was  performed  on  a  MASSCOMP 
MC-5500  scries  minicomputer. 

Cross  Correlation  Analysis  Method 

The  purpose  of  this  analysis  was  to  determine  if  the 
flow  direction  immediately  downstream  of  the  separation 
shock  could  be  deduced  from  wall  pressure  fluctuations. 
Since  the  direction  of  the  undisturbed  boundary  layer  flow 
upstream  of  the  shock  is  known,  the  moving  shock  wave  is 
the  instantaneous  upstream  boundary  of  the  flowfield  where 
flow  direction  near  the  wall  is  not  known.  Only  those  data 
corresponding  to  flow  downstream  of  the  shock  were  needed 
in  order  to  analyze  the  flow  in  this  region.  Therefore,  the 
first  step  was  the  development  of  an  algorithm  to  isolate  that 
fraction  of  the  pressure  signal  corresponding  to  flow 
downstream  of  the  shock.  Two  methods  were  devised  and 
tested,  both  yielding  similar  results. 

In  both  techniques,  a  two-threshold  method  was  used. 
In  all  cases,  the  upstream  transducer  signal  was  used  to 
determine  when  both  transducers  were  downstream  of  the 
shock.  If  the  upstream  transducer  is  downstream  of  the  shock, 
it  follows  that  the  downstream  transducer  must  be  also.  The 
first  threshold,  PT1,  was  "eyeballed"  at  an  estimated  pressure 
value  that  would  represent  typical  pressure  levels  downstream 
of  a  shock.  PT1  was  also  large  enough  to  exclude  the  initial 
pressure  rise  (rising  edge)  due  to  the  shock  passage  (Figure  3). 
The  second  threshold,  PT2,  was  set  at  a  pressure  value  too  low 
to  be  considered  a  value  downstream  of  the  shock.  PT2  was 
also  eyeballed  and  is  also  shown  in  Figure  3. 


Figure  3  Pressure  Time  Histories  and 
Conditional  Algorithm  Thresholds 

Each  shock  passage  is  different  so  the  thresholds, 
particularly  PT1,  for  a  given  run  can  sometimes  result  in  valid 
data  being  discarded.  Or  it  may  start  the  extracted  data  set  too 
soon  and  include  data  with  shock  passage  fluctuations. 
However,  the  results  are  not  particularly  sensitive  to 
physically  reasonable  threshold  settings.  A  brief  discussion  of 
threshold  sensitivity  is  presented  in  the  appendix. 

In  summary,  the  first  algorithm  used  just  the  two 
thresholds  alone  to  find  data  downstream  of  the  shock.  The 
second  method  added  one  further  constraint.  Once  the  level 
of  the  signal  exceeded  the  first  threshold,  each  data  point  was 
examined  to  find  the  maximum  on  the  rising  edge  of  the 
signal.  This  is  referred  to  as  the  "top-finding  "  algorithm.  Once 
past  the  maximum,  the  data  were  assumed  to  be  associated 
with  flow  downstream  of  the  shock  until  values  dropped 
below  PT2,  indicating  the  passage  of  the  shock  downstream. 
This  approach  removes  some  of  the  inflexibility  of  having  two 
"hard-set"  threshold  values. 

The  technique  for  finding  flow  direction  relies  on 
the  cross  correlation  of  two  fluctuating  pressure  signals 
from  two  streamwise  transducers  located  relatively  close  to 
each  other.  The  cross  correlation  equation  is  shown  below: 

Rpp(t)  =  I  p1(t)*p2(t+x) 

where:  p  j(t)  is  the  upstream  channel  of  pressure  data 

P2(t)  is  the  downstream  channel  of  pressure  data 
t  is  the  time  delay  between  channels. 
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In  general,  for  continuously  sampled  data,  cross  correlations 
are  not  calculated  as  shown,  but  are  computed  using  Fast 
Fourier  Transform  (FFT)  algorithms.  FFT  methods  were 
developed  to  reduce  the  computation  time  needed  for  the 
analysis  of  large  data  sets.  Typically  512  or  1024  point 
transforms  are  calculated  and  the  correlation  coefficients  are 
averaged  over  many  records.  The  time  savings  in  using  a  FFT 
algorithm  on  a  small  data  sample  is  not  nearly  as  great. 
Also,  the  algebraic  calculation  is  much  simpler  to  code 
especially  when  the  calculation  data  set  size  needs  to  be 
flexible,  as  was  the  case  here,  as  described  below.  Thus  the 
cross  correlations  in  this  analysis  were  calfculated 
algebraically. 

The  mechanics  of  the  calculation  of  cross  correlations 
downstream  of  the  shock  proceeded  as  follows.  The  number 
of  data  points  (NI)  desired  for  each  "conditionally  extracted 
analysis  data  set”  (CEADS),  typically  32,  64,  or  128  points  is 
input  to  the  code.  The  number  of  points  chosen  depends  on  the 
intermittency  of  the  data.  Data  at  higher  intermittencies  have 
longer  continuous  blocks  of  data  behind  the  shock,  thus 
allowing  longer  CEADS.  Lower  intermittency  values  require 
shorter  CEADS.  A  counter  was  set  at  the  beginning  of  the  data 
file  on  the  upstream  channel  and  marched  through  the  datafile. 
When  the  first  threshold  (PT1)  criterion  was  satisfied,  a 
second  counter  was  started.  NI  data  points  past  the  first 
counter’s  position  were  examined  to  see  that  all  data  points 
satisfied  the  second  threshold  (PT2).  If  a  data  point  violated 
the  second  threshold  (i.e.  P(t)  is  less  than  PT2),  the  first 
counter  was  set  to  the  second  counter's  position,  and 
searching  began  again.  If  no  data  points  violating  PT2  were 
encountered  and  the  second  counter  reached  the  user-set 
value,  the  data  between  the  two  counters  formed  a  new 
CEADS.  This  data  set  was  algebraically  cross  correlated  with 
the  corresponding  data  points  from  the  downstream 
transducer.  Once  the  cross  correlation  calculations  were 
complete,  the  coefficients  were  normalized  by  the  product  of 
the  RMS’s  of  the  CEADS  for  each  channel  in  the  data  set  and 
added  to  a  cumulative  results  array.  The  data  set  counter  was 
advanced  by  one,  the  first  counter  was  moved  to  the  end  of  the 
current  CEADS,  and  the  searching  process  was  started  over. 
When  all  data  were  analyzed,  the  results  array  was  normalized 
by  the  number  of  CEADS  analyzed. 


boundary  layer.  The  single  maximum  at  positive  tq  is 
generated  by  turbulent  eddies  traveling  downstream.  From  tq 
and  the  transducer  spacing  the  broad  band  convection  velocity 
of  the  pressure  carrying  eddies  can  be  calculated.  Its  value  of 
0.67  (496  m/s),  obtained  by  interpolation  of  the  data 

points  bracketing  the  maximum,  agrees  with  previous  work 
[12].  In  contrast,  downstream  of  ’S,’  the  separation  line  from 
the  surface  tracer  experiments,  two  maxima  in  Rpp  are  evident 

and  correspond  to  two  physical  phenomena  (Fig.  5).  The 
Rpp  at  positive  time  delay  corresponds  to  pressure 
fluctuations  due  to  eddies  in  the  separated  shear  layer  flowing 
downstream.  The  broadband  time  delay  for  maximum 
correlation  of  these  structures  is  t=0.006-0.008  msec,  giving  a 
downstream  convection  velocity  in  the  range  365-487  m/sec. 
This  velocity  is  less  than  that  in  the  undisturbed  boundary 
layer  since  the  flow  has  gone  through  the  shock  wave.  The 
second  phenomenon  is  backflow  in  the  recirculating/vortical 
separated  structure.  This  backflow  generates  the  peak  at  t=- 
0.016  to  -0.018  msec,  corresponding  to  a  broadband  upstream 
velocity  of  162-183  m/sec. 
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Figure  4  Standard  Cross  Correlation  of 
Undisturbed  Turbulent  Boundary  Layer 


Discussion  of  Results 

For  reference  purposes  when  discussing  the  cross 
correlations  from  the  conditional  analysis  algorithm  in  the 
intermittent  region,  Figure  4  shows  a  standard  1024  point  FFT 
cross  correlation  result  from  400  records  of  data  taken  with 
two  pressure  transducers  located  in  the  undisturbed  turbulent 
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Figure  5  Standard  Cross  Correlation  of 
Flow  Downstream  of  ’S’ 

As  discussed  earlier,  the  conditional  analysis  space 
time  correlation  code  was  written  to  analyze  signals  in  the 
intermittent  region.  Figure  6  shows  correlations  at  five 
stations  in  the  intermittent  region  for  the  1.9  cm  diameter 
cylinder.  The  legend  indicates  the  normalized  streamwise 
station  of  each  transducer,  the  intermittency  of  the  upstream 
channel  and  the  number  of  CEADS  used  to  calculate  each 
curve.  At  all  values  of  intermittency,  the  "double  peak" 
characteristic  of  separated  flow  can  be  seen.  The  positive  and 
negative  values  of  T  at  which  these  maxima  in  Rpp  occur  are 
the  same  in  all  five  cases.  Based  on  the  transducer  spacing  and 
t ,  the  broadband  upstream  and  downstream  velocities  are  in 
the  ranges  of  195-292  and  292-584  m/sec  respectively.  The 
boundaries  of  these  ranges  are  calculated  using  the  data  points 
bracketing  the  peaks.  From  linear  interpolation  the  actual 
maxima  appear  to  be  about  mid-way  between  these  points, 
giving  upstream  and  downstream  velocities  of  234  and  390 
m/sec  respectively.  Further,  taking  into  account  the  lower 
sampling  rates  for  the  data  sets  in  the  intermittent  region  and 
the  reduced  time  resolution  of  the  cross-correlations,  these 
velocities  are  the  same  as  those  indicated  by  Figure  5.  These 
results  show  clearly  the  flow  structure  downstream  of  the 
instantaneous  shock  location  is  the  same  at  all  stations  within 
the  intermittent  region.  Further,  it  is  the  same  as  the  flow 
structure  in  the  separated  flow  downstream  of ’S.’ 

Figure  7  shows  cross  correlations  using  the  1.27  cm 
diameter  cylinder.  All  three  stations  show  the  characteristic 
"double  peak",  indicating  backflow  immediately  behind  the 
shock.  No  variation  in  time  delay 
occurs  at  t*-0.014  msec  and 
upstream  and  downstream  velocities  of  209  and  292-365 


at  R  is  seen.  R 
t=0.08^6l  msec  givTng 


m/sec  respectively.  Examination  of  Figures  6  and  7  reveals 
that  both  cylinders  have  the  same  time  delays  for  both  peaks, 
within  the  sampling  rate  accuracy,  indicative  of  the  same  flow 
structure  immediately  downstream  of  the  shock. 
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Figure  6  Cross  Correlations  of  CEADS- 1.9  cm  Cyl. 


Figure  7  Cross  Correlations  of  CEADS- 1.27  cm  Cyl. 
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Conclusions 


A  method  of  detecting  flow  direction  downstream  of 
the  unsteady  separation  shock  using  fluctuating  wall  pressure 
signals  has  been  developed  and  tested  in  a  Mach  5  shock 
wave  turbulent  boundary  layer  interaction  induced  by  circular 
cylinders.  In  this  method,  flow  direction  is  deduced  from  cross 
correlations  of  that  part  of  the  pressure  signal  downstream  of 
the  moving  shock  wave.  The  results  show: 

(1)  For  all  cylinder  flows  examined  the  cross  correlations 

A. 

show  a  maximum  Rpp  at  negative  time  delay  indicating  that 
backflow  exists  immediately  downstream  of  the  shock  at  all 
stations  in  the  intermittent  region. 

(2)  Separation  essentially  occurs  across  the  shock  in  the 
intermittent  region  upstream  of  unswept  circular  cylinders  and 
hence  the  separation  point  itself  undergoes  a  large  scale 
streamwise  motion.  It  appears  thaf  the  well  defined  separation 
line  from  the  kerosene  lampblack  pattern  delineates  the 
downstream  boundary  of  a  region  of  intermittent  separation. 
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Appendix 

Conditional  sampling/analysis  algorithms  are  a 
relatively  new  type  of  analysis  being  used  in  these  typies  of 
interactive  flows.  When  data  is  converted  into  a  new  signal, 
box  cars  (0’s  or  l’s)  for  shock  motion  analysis  for  instance,  or 
split  into  different  segments  based  on  given  criteria,  the 
sensitivity  of  the  method  to  the  decision  criteria  should  be 
investigated.  This  was  recently  discussed  by  Dolling  and 
Brusniak[6]  regarding  box  car  analysis.  To  confirm  the 
validity  of  the  conditional  analysis  space-time  correlation 
code,  a  threshold  sensitivity  analysis  was  performed. 

Figure  8  shows  the  conditional  analysis  cross- 
correlation  results  for  different  thresholds  for  the  same  set  of 
data.  Two  cases  are  shown  for  clarity.  Although  the 
magnitudes  of  the  coefficients  have  changed,  the  positive  and 
negative  values  of  R  remain  at  the  same  time  delays  (1). 

**  mix 

This  was  also  observed  for  several  other  threshold 
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combinations.  Thus,  although  the  values  of  R  fluctuate, 

PPn 


the  corresponding  values  of  I  are  insensitive  to  changes  in 
threshold  values. 
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Figure  8  Sensitivity  of  Cross  Correlations  of  CEADS 

A  similar  result  is  obtained  when  comparing 
algorithms.  Figure  8  shows  results  for  both  algorithms  tested. 
Each  case  had  the  same  thresholds.  The  top  finding  algorithm 
had  fewer  CEADS  which  met  the  32  point  analysis  data  set 
requirement  due  to  the  stricter  requirement  of  starting  the  data 
set  with  the  maximum  point  past  the  first  threshold,  PT1.  As  a 
result,  more  fluctuations  due  to  the  shock  passage  were 
eliminated,  and  the  corresponding  correlation  values  due  to 
flow  behind  the  shock  are  slightly  higher.  However,  the  time 
delay  values  for  the  peaks  remain  unaffected. 
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Appendix  B 


Flowfield  Model 

The  reason  why  the  surface  *streak  lines  are  in  the  downstream  direction 
upstream  of  ’S’  can  be  explained  using  a  relatively  simple  model.  Consider 

the  flow  on  centerline.  If  the  instantaneous  surface  shear  stress  in  the 

/ 

intermittent  region  is  modelled  as  a  step  function  (Fig.  13),  then  the  mean 
wall  shear,  x„,  which  the  surface  streaks  respond  to,  is  given  by 
Xh-=(1-7)  Xu  +YX d 

where  xu  and  xd  are  the  average  wall  shear  stresses  in  the  upstream  and 
downstream  zones.  xw  will  equal  zero  when 

T-  1  =  Jd_ 

'i  V 

Since  T„  is  the  wall  shear  stress  of  the  incoming  supersonic  boundary  layer 
and  x d  is  at  the  upstream  boundary  of  the  separated  flow,  the  ratio  xu/x d, 
given  by  N  will  be  large.  Hence,  xw  =  0  when  y  =  N/(l+N).  For  large  values 
of  N,  y  is  close  to  1.  Thus  the  mean  shear  stress  at  the  wall  can  be  in  the 
downstream  direction  even  when  the  flow  is  separated  for  the  major  fraction 


of  the  time. 
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Figure  13  -  Shear  Stress  Model 


